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Abstract. Since the first observations by Kaufmann et al. (1970), special attention has been paid to static pressure-balanced 
structures in the form of magnetic holes or humps observed in regions of the solar wind and of planetary magnetosheaths 
where the /3 parameter is relatively large and the ion perpendicular temperature exceeds the parallel one. Although alternative 
interpretations have been proposed, these structures are usually viewed as associated with the mirror instability discovered in 
1957 by Vedenov and Sagdeev. After reviewing observational results provided by satellite missions, high-resolution numerical 
simulations of the Vlasov-Maxwell equations together with asymptotic and phenomenological models of the nonlinear 
dynamics near the instability threshold are discussed. The constraining effect of the mirror instability on the temperature 
anisotropy associated with a dominant perpendicular ion heating observed in the solar wind is reported, and recent simulations 
of this phenomenon based on an elaborated fluid model including low-frequency kinetic effects are briefly mentioned. 
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1. INTRODUCTION 

This paper is based on three one-hour lectures devoted to the mirror structures observed in space plasmas. Being asso- 
ciated with the nonlinear development of the mirror instability, these structures are often referred to as nonlinear mirror 
modes. Their understanding results from situ measurements performed by satellite missions, direct numerical simula- 
tions of the Vlasov-Maxwell equations, and also theoretical developments based on asymptotic and phenomenological 
modeling. An overview of the main conclusions of these approaches is presented, including a discussion of enlighten- 
ing figures published in the literature together with an extended bibliography. 

For various reasons, mirror modes deserve a special interest. They are associated with a fundamental instability 
in collisionless plasmas, whose nonlinear saturation is still not fully understood. They also benefit from accurate 
observational data originating from various space plasma environments, and are expected to have major effects on the 
global dynamics of the plasma both at the level of the transport properties and as constraining temperature anisotropies. 

Conspicuous examples of mirror structures in the form of magnetic enhancements (usually referred to as peaks or 
humps) or depressions (holes or dips) are provided by Figs. 1 of references [1] and [2] that report measurements in the 
Earth's and the Jovian magnetosheaths respectively. Such structures are in fact ubiquitous in spaces plasmas such as the 
solar wind and the planet magnetosheaths, in regions where the perpendicular temperature of the ions dominates the 
parallel one. The mechanism leading to their formation is however not fully understood and could also be non unique. 
It should furthermore be noted that the same temperature anisotropy can also lead to the ion cyclotron instability that 
displays significantly different properties ([3], Section 7.2). 

The question also arises of the origin of this temperature anisotropy. In the planetary magnetosheaths [4], large 
perpendicular to parallel ion temperature anisotropy ratios are expected in two regions (see Fig. 1 of JU), namely 
downstream of the quasi-perpendicular bow-shock (i.e., the region where the angle between the interplanetary mag- 
netic field and the shock normal exceeds 45°) where a preferential perpendicular heating of the plasma is due to shock 
compression, and also close to the magnetopause where field-line draping is important, the anisotropy being caused 
by flow of plasma with high parallel velocity out of the ends of the magnetic flux tubes. In the solar wind, the problem 
is not really settled |@] and will be briefly addressed in the following. 

An additional difficulty in understanding the formation of mirror structures originates from their presence in 
regions where the plasma is linearly stable. As discussed below, a possible interpretation relates this observation 
to the subcritical character of the bifurcation. Alternative interpretations have nevertheless been suggested. It was in 
particular claimed that mirror structures are in no way related to the mirror instability but should rather be viewed as 
slow-mode magnetosonic soli tons [7]. Differently, the mirror instability was viewed as the driver for the developments 
of magnetic dips and humps that evolve differently when convected away from the region where the instability is 



triggered H], making the observed structures remnants of mirror modes created upstream of observation point fljj. 



2. THE MIRROR INSTABILITY 



The mirror instability was discovered in 1957 by Vedenov and Sagdeev |10y who performed, in the case of a 
proton-electron plasma, a linear perturbation analysis of the Vlasov-Maxwell equations about a bi-Maxwellian ion 
distribution function under the assumption of cold electrons (a regime for which the parallel electric field is negligible), 
leading to the instability condition 
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where Bq is the ambient field, no the plasma density, Tt and 7ji the ion perpendicular and parallel temperatures. The 
same instability condition is given in the appendix of 111 1|] w hich deals with a cylindrical plasma in the so-called pinch 
configuration. Interestingly, the volume containing Ref. [ 10] also includes an analysis by Rudakov and Sagdeev fl2l1 of 
an equivalent instability based on a bi-fluid description, which leads to a significantly different threshold. The name of 
the instability was suggested by the magnetic bottles it creates (Fig. 3.9 of IU3IP . and by the role of the mirror force in 
the instability development. Quoting 01211 . "let some perturbation of the uniform density be applied. Since the plasma 
is diamagnetic, the field is reduced where the density is increased. In a non-uniform field, however, a force — /xgradB 
acts on the plasma, and this increases the original perturbation." A more quantitative phenomenology discriminating 
between resonant and non resonant particles is presented in IU4I1 . 

Mirror instability has an essentially kinetic origin, in spite of the existence of a corresponding magnetohydrody- 
namic instability. When the equilibrium state is homogeneous, it is a zero-frequency instability in the rest frame of 
the plasma. As a consequence, the mirror modes are not propagating in this frame. Differently, the so-called drift mir- 
ror instability that develops in the presence of mean-field and or density gradients is oscillatory and involves a finite 
propagation velocity ifTBl Il6i [nil . 

The growth rate of the mirror instability was first evaluated in the large-scale (or quasi-hydrodynamic) limit lfl8ll . A 
main observation is that the instability is driven by the Landau wave-particle resonance, which explains the discrepancy 
with the bi-fluid approach. It is quenched at small scales by finite Larmor radius (FLR) effects [ 1 9, 20] . Near threshold, 
for an electron-proton plasma, under the assumption of cold electrons and of a bi-Maxwellian distribution function for 
the ion velocities, the growth rate is given by 
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where j3|| = fi±T^/Tj_ is the parallel ion beta, v th || = (7j| jm) l l 2 and Vthj_ = {T^/m) 1 ! 2 are the parallel and perpendicular 
thermal velocities of the ions, and = v^/Q. = (1 /Q.)(Tj_/m) 1 / 2 denotes the ion Larmor radius. Here Q = eBo/(mc) 
is the ion gyrofrequency, e and m being the charge and the mass of the proton, and c the speed of light. Note that 
the factor \k^ | (with the absolute value) reflects the Landau resonance. The growth rate expression is exemplified 
for different temperature anisotropy ratios and different beta parameters in Fig. 1 of H20fl . Equation Q shows that 
the mirror instability dominantly affects wavevectors substantially oblique to the ambient field, and can thus be 
viewed as quasi-transverse near threshold. Furthermore, it is associated with magnetic fluctuations predominantly in 
the direction of the background field (quasi-parallel linearly polarized modes). The difference is conspicuous with 
the ion cyclotron instability that develops under the same temperature anisotropy conditions but that has a non- 
zero frequency, a maximum growth rate in the direction parallel to the ambient field and predominantly amplifies 
perpendicular magnetic-field components. 

When the electrons are warm with a temperature anisotropy, the parallel electric field is relevant and the instability 
condition becomes JHHHH (see also Eq. (49) in Section 9.8 of 121 ). 
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Plots of the growth rate versus the wave number are presented in Fig. 2 of [25] and Fig. 10 of j26ll . 

The mirror instability in the case of cold electrons and a general distribution function for the ions is considered in 
Il8ll27il28il29il23ll . For a distribution function of the form /(vjj, VjJ, the instability condition then reads (ps = Bq/Stc) 
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We also mention a gyrokinetic description of the mirror instability presented in [30]. 



3. SATELLITE OBSERVATIONS 

Mirror modes are observed in various space-plasma environments displaying an ion temperature anisotropy of the 
form Tj± > 7}|| and large or moderate j3 parameters, including regions where the plasma is linearly stable. They have 
been reported in the Earth's magnetosheath 13 ll l32l 33 l 34 l l35l l36l 111 1371 13811 and the downside magnetosphere 13 911 . 
the Jovian Holllllll], Saturn's BHU], and Venus'J|44 L H magnetosheafhs, the Io wake JH,li3], the comets Halley 
lf48h and Giacobian-Zinner ll49ll . the solar wind |^jl50l 5 ill, the heliosheath 11521 l53ll . the interplanetary space lf54Tl . or 
ahead of interplanetary coronal mass ejection B55I1 . Such a frequent occurrence of mirror structures, their relation with 
a basic plasma instability and also their potential influence on the global dynamical processes such as the diffusion of 
energetic particles 115611 . the transport properties of high j3 astrophysical plasmas 15711 or the constraint on temperature 
anisotropy (see Section 7), motivate a detailed analysis of their formation and time evolution. 

Mirror structures are characterized by a few important properties. They are quasi-static in the plasma frame [58] 
when the plasma is homogeneous (in contrast with MHD modes such as slow modes). Drift mirror modes (that are 
propagating) indeed require density or ambient-field gradients. Mirror modes involve only a small change in the 
direction of the magnetic field, being dominated by the variation of its amplitude. They are thus referred to as linear 
structures. The magnetic field fluctuations mostly affect the parallel component and are thus almost linearly polarized 
in the direction of the ambient field. The relation between linear polarization and the non-propagation of the mirror 
modes is discussed in ll59ll . The spatial variations mostly take place in highly oblique directions (typically 70° or 80°) 
and are thus qualified as quasi-transverse. Furthermore, Cluster multi-spacecraft observations indicate that mirror 
modes can be viewed as cigar-like structures, quasi-parallel to the ambient field with a transverse scale of a few ion 
Larmor radii ||38[ l60ll . A important property is the anticorrelation between the variations of the density and of the 
magnetic field amplitude. It is exemplified in the case of nearly sinusoidal modes and of magnetic humps observed in 
the Earth's magnetosheath, respectively in Figs. 1 and 2 of 16111 . The anticorrelation in the case of magnetic holes is 
reported in Fig. 3 of [62]. This figure also displays the fluctuations of the plasma and the magnetic pressures, together 
with their sum that turns out to be almost constant, thus indicating a pressure-equilibrium regime. 

Magnetic depressions with various shapes are observed in different space-plasma environments (Fig. 1 of [63], Fig. 1 
of Il64ll or Fig. 2 of 15011 ). and the question arises of the possibly different origins of these structures. Furthermore, huge 
pressure-balanced structures in the form of magnetic holes with cross-sectional width of hundreds or even thousands 
proton gyroradii were reported in the solar wind at 1 AU (Fig. 3(b) of [65]). Temperature anisotropy and magnetic 
field measurements give direct evidence that the observed magnetic holes are stable remnants of magnetic pressure 
depletions generated in a source region closer to the Sun, likely through the mirror instability, and subject to the solar 
wind expansion. It was also suggested on the basis of hybrid numerical simulations that in the high speed solar wind 
magnetic holes can be generated by large-amplitude right-hand polarized Alfvenic wave packets propagating at large 
angle to the ambient magnetic field [66] . 

Various sophisticated methods (discussed in 1 60] ) were developed to identify observed low-frequency modes ll67ll 
and among them mirror modes. Indeed, slow and mirror modes both display antic orrelated magnetic and density 
fluctuations and mirror and ion-cyclotron modes grow under the same temperature anisotropy condition. Nevertheless, 
in the Earth's magnetosheath, heavy ions (He +2 , also called a lpha pa rticles) tend to make the mirror instability 
dominant for j3 > 1 and the ion cyclotron instability at smaller j3 1681 |25i |69l H F70I p7ll I7I1 . 

Mirror structures are observed in the form of either magnetic depressions (dips or holes) or enhancements (peaks 
or humps), as already conspicuous in the first mirror mode observations reported in Fig. 1 of 13111 . Figure 2 of 14311 
shows, in the case of the Saturn's magnetosheath, that magnetic peaks are mostly observed near the bow shock while 
magnetic holes are present near the magnetopause. A consistent observation was made in the Jovian magnetosheath 
J2] where "mirror structures are present everywhere [•••] with peaks appearing in the middle of the sheath and holes 
close to the magnetopause, at least for the day side magnetosheath which represents most of the data." Similarly in the 



Earth's magnetosheath i37ll . "peaks are observed closer to the shock. Holes may be observed everywhere in the sheath 
with a preference for the vicinity of the magnetopause and for mirror stable conditions." 

The question then arises of the conditions prescribing the shape of the mirror structures. A relation with the 
amplitude of the ambient magnetic field in the Earth's magnetosheath was first mentioned in ll34ll . Figure 2 of this 
article shows magnetic peaks encountered in low-field regions, dips in high-field regions and also near-sinusoidal 
waveforms. More specifically, Fig. 8 of [2], that displayed the variations of the magnetic field amplitude and of the j3 
parameter in the Jovian magnetosheath, suggests that magnetic holes are favored by low values of j3 and humps by high 
values. Quantitative information is provided by Fig. 6 of the same reference that displays the occurrence distribution 
of peaks and holes versus j3. A strong concentration of dips is observed for j3 in the range (0.5, 1), that decays to 
zero by /3 ss 3.5. Peaks are mostly observed between j3 = 2.5 and 6, with a maximum near 4.5. The distribution of 
quasi-sinusoidal modes is rather flat, between /3 = 0.5 and 5. A quantitative characterization of the shape of the mirror 
structures is provided by the skewness of the magnetic fluctuations IB7I1 or by their peakness OJ] (which is in practice an 
essentially equivalent concept defined as the skewness of the time series representing the total wavelet content between 
two chosen scales of the original magnetic field fluctuations). Typically, a positive skewness or peakness corresponds 
to peaks, and a negative one to holes. Figure 3 of 13711 and Fig. 2 (left) of [1] show the skewness or the peakness 
versus j3ii for all the mirror structures detected in a given period of time in the respective observations. In spite of the 
dispersion of the clouds of points in both pictures, a correlation between the value of J3| and the sign of the skewness 
or the peakness is conspicuous: larger j3ii favors magnetic humps and lower j3ii magnetic holes. 

In addition to the value of J3||, another parameter that can affect the mirror mode features is the ion temperature 
anisotropy ratio Zj_/Tj| . Figure 3 of [ 1] shows, in the plane (J3|| , 7j_/7j|), the distribution of mirror modes of different 
types (specifically, those with a peakness smaller than —0.6 and those with a peakness larger than 0.3). A main 
conclusion is that "peaks are typically observed in an unstable plasma, while mirror structures observed deep within 
the stable region appear almost exclusively as dips." Consistent with this statement, Fig. 6 of [65] displays, in the same 
parameter plane, the distributions of magnetic holes observed in the solar wind at a distance of 1 AU. These structures 
are almost all located in the stable region, deeper holes appearing to congregate near the instability threshold. Indeed, 
as noted in [73], the solar wind is almost stable against the mirror instability, leading to the authors to conclude that 
Ulysses satellite "observed structures generated by the mirror mode instability, which remained after the distribution 
relaxes to a marginally steady state." So, "although the plasma surrounding the holes was generally stable against 
the mirror instability, there are indications that the holes may have been remnant of mirror mode structures created 
upstream of the point of observation." 

A detailed analysis of the relation between the shape of the mirror structures and the distance to threshold, estimated 

by the mirror parameter Cm = j8,j_ — lj (a description strictly valid for cold electrons and bi-Maxwellian ions but 

viewed as satisfactory, given the accuracy of the data), is presented in Ref. i37ll . Figure 6 shows that magnetic holes 
(and negative skewness) mostly occur when the plasma is linearly stable (Cm < 1). Figure 7 shows that magnetic peaks 
(and positive skewness) are observed in a regime well above threshold (Cm > 1)- A comprehensive overview is given 
by Fig. 5 that displays the average skewness as a function of Cm- Below threshold the skewness is negative, while a 
positive skewness is obtained only above threshold. The figure also includes an insert that indicates the tendency of the 
Earth's magnetosheath plasma to be mostly observed in a marginally stable state with respect to the mirror instability, 
in agreement with other publications. It furthermore turns out that peaks are a minority among the mirror structures. 
In yj], it is reported that 14% of the observed structures are magnetic peaks and 19 % are holes. Similarly in 
18.7% of the observed mirror structures are peaks and 39.7% are holes. In agreement with QJJ] and B43I1 . observations 
reported in [37] indicate that holes are preferentially observed close to the magnetopause and peaks are more frequent 
in the middle of the magnetosheath. In the solar wind also, magnetic holes occur much more frequently than magnetic 
enhancements B74I1 . 

4. A MHD-LIKE MODEL FOR STEADY MIRROR STRUCTURES 

Although the mirror instability is driven by kinetic effects, some properties of stationary mirror structures can be 
captured within the anisotropic magnetohydrodynamics (MHD) supplemented by suitable equations of state. The 
simplest anisotropic MHD model is provided by the double adiabatic or CGL approximation, for Chew Goldenberg 
and Low ([75], see also B76I1 ). It consists in closing the fluid hierarchy derived from the Vlasov equation by assuming 
no heat fluxes and concentrating on scale large enough for the Hall effect to be negligible, which leads to the equations 
of state px = pB and pii = p 3 /B 2 , where B is the magnitude of the magnetic field and p the plasma density. This 



approximation is valid in regimes where typical phase velocities are much larger than the particle thermal velocities, but 
is not appropriate for mirror modes that are static in the plasma frame. The big discrepancy between the kinetic theory 
and the CGL prediction, which concerns not only the mirror growth rate but also the instability threshold (namely of 
a factor 6), is discussed in IT76ll and numerically exemplified in Fig. 5 of lf77ll . Ad hoc anisotropic polytropic equations 
of state with parallel and perpendicular indices chosen in such a way that the actual mirror instability threshold is 
reproduced, have also been suggested ( i78l,l79ll and references therein). 

An isothermal or a static limit is in fact more suitable for mirror modes. A series of equations can be derived for 
the gyrotropic components of the even moments, and using the assumption of bi-Maxwellian distributions, simple 
equations of state can be obtained, which predict the correct threshold of the mirror instability. This model B80I1 is 
potentially relevant for describing the large-scale features of steady mirror structures, since Landau damping vanishes 
for static solutions and FLR-corrections are negligible at large scales. These effects are however needed in order to 
reproduce the correct instability growth rate and to capture the time dynamics (see Section 7.3). 

To construct the model of Ref. ifsotl ^ one assumes a static regime characterized by a zero hydrodynamic velocity and 
no time dependency of the other moments. Projecting the ion velocity equation along the local magnetic field (whose 
direction is defined by a unit vector b) leads to the parallel equilibrium condition 

b m ^-P mn = (6) 
ox„ 

for the (gyrotropic) pressure tensor Pjj — p±rtjj + p\\Tjj where rey = 5,-y — bjbj and Ty = bfij are the fundamental 
gyrotropic tensors. 

Proceeding in a similar way at the level of the equation for the heat flux tensor after contraction with the two 
fundamental tensors ny and Ty, gives the conditions 

f d 

b, nij —R ijkl =0 (7) 

brtj-^Rij H =0, (8) 

where the (gyrotropic) fourth-order moment takes the form 

Rijki = r\\\\TijTki + r\\ ± y{nijZ kl } + r xx y{n i} T k i}. (9) 

Here 5^ indicates symmetrization with respect to all the indices. The closure then consists in assuming that the fourth 
order moments are related to the second order ones as in the case of a bi-Maxwellian distribution, namely ry | = 3p| /p , 

1|± = P\\Pi-/P an d r ±± = 2p]_/P- Straightforward algebra then leads to the following equations of state for the 
temperatures normalized by their equilibrium values 
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where A = (r) 0) /^ (0) ) - 1 measures the anisotropy of the equilibrium state. Equation ( fTTT ) also rewrites 

7j_ = 7j| f(B/Bo) with f(u) = u/[u— A(A + l) -1 ]. Similar equations of state were derived in i8lll . using a 
fully kinetic argument. 

With the above equations of state, the problem of stationary structures is amenable to a variational formulation. 
Indeed, estimating the mechanical work produced by a small expanding plasma cylinder and using the conservation of 
the particle number n and of the magnetic flux, one computes the effective energy per particle as E = 7j| [F{B/Bq) +lnn] 

with F(B/B ) = J^ B ° ^r^ du = H( A + l)( B / B o) - A] - ]n(B/B ). At equilibrium, the total (magnetic + elastic) 
potential energy of the system — j[(B 2 /%%) +nE(n,B)]dr should be minimum, taking into account the constraints 
of particle conservation and the Lagrangian nature of the frozen-in magnetic flux. This leads one to consider the 
functional 
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where J3|| = S7Cp^' /Bq. Minimizing this quantity under the above constraints provides the equation determining the 
possible stationary plasma configurations in the form 
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which expresses that zero force is acting on each fluid element. It in particular follows thatB • V8 Jf? /8n = 0, indicating 
that 8J€ '/ 8n is a constant (denoted by In C) along each magnetic flux line. In the special case where this constant is 
the same for all the lines, one gets in the case of cold electrons the relation n = Cexp[—F(B/Bo)]. 

In two space dimensions, the magnetic field is written as B = (—YyiYx) m terms of a stream function y/(x,y). 
Substituting in Eq. Sl3[ and linearizing about a uniform background by writing C = 1 and Vi// = e + Vv where e is the 
unit vector along the x-axis and v(x,y) a small quantity, one gets the equation 

[2/j3|| -2A(A + l))]v„+ [2//3|| +A]v, y = 0. (14) 

This equation is hyperbolic provided A (A + 1) — 1/j8ii > 0, or equivalently A — 1 > 1/j8j_ > 0, which reproduces the 
instability condition (fl} for mirror modes. 

When concentrating on one-dimensional solutions by assuming that the fields depend on a single Cartesian coor- 
dinate x in the formB = (cos 0, sin 9 + <j)(x),Q) and n — n(x), Eq. ([T3T l has only piecewise constant solutions (with 
specific jump conditions). When assuming periodic boundary conditions, they take the form 

Bi = (cos0,sin0-a) m = l + 5/Z> if 0<x<X (15) 
B 2 = (cos 6, sin0 —b) n%=\ — 8ja if A < jc < 1 , (16) 

that displays a discontinuity at x — (or 1) and at x = A = b/(a + b). The expression of A, n\ and ni are prescribed by 
the constraints that the mean magnetic field and the number of particles are equal to those of the unperturbed solution. 
Density humps (and magnetic holes) correspond to A < 1/2, while density holes (and magnetic humps) are associated 
with A > 1/2. The positive parameters 8, a and b are then determined in terms of j8n, the temperature anisotropy 
and the angle 9, by minimizing the energy of the solution. It turns out that in addition to the trivial stationary point 
a = b = 0, the system can admit non trivial solutions associated to mirror structures. Details are presented in [80] where 
Figs. 4 and 5 display profiles of the magnetic field and of the density for 6 = 70°, 7j_/Tj| = 1.25 in the cases = 2 
and J3|| = 16 respectively. In the former case which corresponds to a subcritical regime, magnetic holes and density 
humps are obtained, while in the later case where the plasma is mirror unstable, the structures are magnetic humps and 
density holes. The variation of the parameter A that characterizes the shape of the structures in terms of the distance 
to the instability threshold for various angles 0, as predicted by this model, is plotted in Fig. 8 of l37ll . Especially 
for 6 = 70°, a noticeable qualitatively good agreement is obtained with Fig. 5 of the same reference that displays the 
skewness of the magnetic fluctuations associated with the mirror structures observed by the Cluster mission, versus 
the distance to threshold. 



5. VLASOV-MAXWELL NUMERICAL SIMULATIONS 

Observational data have the great advantage of referring to the "real world." But space plasmas are complex systems 
involving a large number of simultaneous phenomena that possibly affect one another. Furthermore, observations can 
hardly provide information on the dynamics. In contrast, numerical simulations can focus on a specific phenomenon 
and permit the detailed study of its time evolution. Thus, in spite of the limitations of the computational facilities, 
numerical simulations of the Vlasov-Maxwell system are most useful for understanding the nonlinear development 
of the mirror instability and the formation of the mirror structures. The present-day computers enable such a study, at 
least in one space dimension, i.e., when the spatial variations are considered only along a direction making a prescribed 
angle 9 (usually associated with the largest linear growth rate) with the ambient magnetic field (taken in the z direction 
in the following). Recent such simulations for an electron-proton plasma are presented in [82], assuming periodic 
boundary conditions in the space variable. Two computational algorithms are used for the kinetic description of the 
ions, the electrons being viewed as a charge-neutralizing fluid (with or without mass) keeping a constant temperature 
(taken equal to zero or very small), and the electromagnetic field being advanced by solving the Maxwell equations. 
Such codes are usually referred to as hybrid. One of the description is the so-called particle in cell (PIC) algorithm 



that essentially consists in solving the Vlasov equation by the characteristic method. It is especially efficient, and is 
thus conveniently used for simulations in large computational domains. A low statistical noise, especially needed for 
simulations near the instability threshold, nevertheless requires a large number of particles per cell. The other method 
consists in an Eulerian description for the Vlasov equation. The advantage of the latter approach is the absence of 
statistical noise, but it requires more computational power than the PIC algorithm, especially near threshold where it 
is necessary to accurately resolve for the velocity variables in order to minimize the numerical diffusion. 

A hybrid PIC simulation of the nonlinear development of the mirror instability near threshold in an extended domain 
(2048 ion inertial lengths), is presented in B82I1 . using a resolution of 1024 cells with 500 000 particles in each cell. 
The electrons are assumed massless and almost cold. The initial ion distribution function is taken bi-Maxwellian with 
J3|| = 1 and 7±/7|| = 1.857, thus corresponding to an initial distance to threshold Fq = j3j_(7j_/7j| — 1) — 1 = 0-6. 
The retained direction is that of the most unstable mirror mode (8 = 72.8°, for which the linear growth rate is 
7 = 5 x 10 _3 £2 P where £2 p is the ion gyrofrequency). Figures 1 and 2 show that the saturation of the linear instability 
leads to the formation of magnetic humps that, as time elapses, undergo a coarsening process which leads to their 
strengthening. As a consequence of the large extension of the computational domain, a large number of Fourier 
modes are excited in this simulation. In agreement with the evolution of the structures, the skewness of the magnetic 
fluctuations rapidly becomes positive, increasing up to the moment when the energy of the magnetic fluctuations 
saturates. An interesting issue concerns the instantaneous distance T to threshold defined by the left hand side of 
Eq. ©, that rapidly departs from the bi-Maxwellian expression, indicating a significant evolution of the functional 
form of the distribution function. It also does not relax to zero but reaches negative values, while the energy and the 
skewness of the fluctuations are still growing under the effect of the nonlinearities. 

The departure from a bi-Maxwellian profile for the ion distribution function, that develops during the simulation, 
is exemplified in Fig. 4 which shows that the space-averaged instantaneous distribution function (/) displays a clear 
flattening near V|| = 0. This effect can be associated with a diffusion process in the velocity space due to quasi-linear 
effects (see Section 6.4). 

It is noticeable that when the same simulation is performed in a much smaller computational domain, quasi-linear 
effects are not significant. In contrast, one observes an oscillation of the energy of the magnetic fluctuations and also 
of the skewness with a frequency which is consistent with the bounce frequency ffl^ = (1 /2)v^ ± k\(8B/B f of the 
trapped particles, suggesting that particle trapping is at the origin of these oscillations (which are not present when the 
computational domain is large enough). 

Similar PIC simulations were previously performed at larger distance from threshold (for which a lower resolution 
is sufficient) [8]. They also clearly indicate that an initial noise in a mirror unstable regime leads to the formation of 
magnetic humps whose number decreases as time elapses (Fig. 4). In addition, these simulations show that in a plasma 
with isotropic ion and electron temperatures (thus mirror stable), a localized magnetic perturbation in the form of a 
finite-amplitude hump relaxes (Fig. 3), while a perturbation in the form of a finite-amplitude magnetic hole persists 
(Fig. 1). 

As already mentioned, Eulerian simulations are also reported in 1 82] in the case of a relatively small domain (307C 
ion inertial lengths). In the supercritical regime (cold electrons and initial bi-Maxwellian ion distribution characterized 
by J3|| = 6 and T±/T^ = 1.25, together with an angle 9 = 83.86°), the nonlinear development of the mirror instability 
leads to the formation of a magnetic hump anticorrelated with a density hole (Fig. 7). A similar evolution is shown in 
Fig. 8 at larger distance from threshold (j8|| = 15, T±/Tn = 1.4 and 9 = 78.53°). Oscillations of the magnetic-energy 
fluctuations are visible (as in the PIC simulations in small domains), that are believed to be the signature of particle 
trapping. 

Below threshold (isotropic ion temperatures, J3ii = 6 with 9 = 83.82°), an initial magnetic hole with no density 
perturbation (Fig. 15) evolves slightly but persists, while an anticorrelated density hump develops. Remarkably, a 
similar evolution takes place above threshold, as exemplified in Fig. 13 with j8|| = 6, 7±/7|| = 1-36 and 9 — 83.82°. It 
is of interest to note that, above threshold, the magnetic profile develops overshoots at the boundary of the depression, 
as in space plasma observations. 

In the above simulations, initial magnetic holes are shown to be persistent, but their formation starting from an 
initial noise is a different issue. Two mechanisms were reported for the development of such structures. In [82], a PIC 
simulation is presented (Fig. 16) in a domain extending on 1024 ion inertial lengths, in conditions relatively far above 
the mirror instability threshold (fin — 1, = 4 with 9 = 50.5°, leading to a linear growth rate of 0.156£2 p ). At 

early time, magnetic humps are formed as in the previously described simulations. However, later on, these structures 
progressively evolve towards magnetic holes. The instantaneous distance T to threshold decays but now remains 
slightly positive. The system is thus continuously stirred and the coarsening is less efficient. In particular, there is 



no formation of isolated structures. Furthermore, the energy of the magnetic fluctuations does not grow monotonically 
up to its saturated value but reaches a sharp maximum before decaying to an equilibrium value (smaller than the 
maximum by more than a factor 2). The skewness has a similar evolution with the same time scale. It reaches a 
(positive) maximum at the same moment as the energy, and then decays to a negative value. It is noticeable that in 
this simulation, the ion distribution function remains relatively close to a bi-Maxwellian with no significant flattening 
near the zero parallel velocity. It turns out that the transition of the magnetic structures from humps to holes does not 
take place with a larger j8|| (e.g., j8|| = 2). Another scenario for hole formation is reported in [37] and concerns a PIC 
simulation in an expanding domain modeling the magnetosheath 18311 . Plots of the skewness versus j8|| and versus the 
(bi-Maxwellian) distance to threshold F* = j3i(7j_/7j| — 1) — 1 as the system evolves, are presented in Fig. 9 of l37ll . 
Magnetic humps evolve into holes as they are convected closer to the magnetosphere within the plasma depression 
layer where the plasma expends and its properties change drastically. In particular jS|| and F* decrease, the transition 
from humps to holes occurring for F* m 1 .7. 

From the combined views of observations and models, the following scenario is proposed in [37] for the evolution 
of the mirror structures in space plasmas. "Large peak structures grow out from moderately unstable plasma (typically 
observed behind the bow shock), reaching saturation at sufficiently large distance from threshold (in the middle 
magnetosheath depending on the convection time), and, as j3 is further decreased closer to the magnetopause, the 
plasma turns to be mirror stable, which is accompanied by a decay of the peaks to the profit of hole structures which 
can survive these conditions." 



6. TOWARD A NONLINEAR THEORY NEAR THRESHOLD 
6.1. A reductive perturbative expansion 

A natural way for extending the linear theory to the weakly nonlinear regime that develops near threshold is to 
perform a long-wavelength reductive perturbative expansion, based on the observation that near threshold the unstable 
modes are localized at large scale [84, 82(1 . 

From the Vlasov equation, one derives the equation for the proton velocity 

d« 1 „ e , 1 

— + -V p (E + -uxB =0, (17) 

at p trip c 

where, for cold and massless electrons, E = — (u — — ") x B, with ;' = (c/47r)V x B. The ion pressure tensor is 

c V neJ 

rewritten as the sum of gyrotropic and non gyrotropic (often called gyroviscous) contributions p = (p±n + p^z) +11, 

with n = I b<E)b and T = b<E) b, where b =B/\B\ is the unit vector along the local magnetic field. Equation ( fTTI i is 
then rewritten in the form 

Only the projection on the plane perpendicular to the ambient field, obtained by contracting Eq. dT8b with the tensor n, 
is actually needed. In order to address the asymptotic regime, assuming the ambient field in the z direction, we rescale 
the independent variables as prescribed by the linear growth rate, in the form X — s/ex, Y = s/Ey, Z = ez, T = S 2 t, 
where e measures the distance to threshold. We furthermore expand any field (p in the form 

<p = £ e"/ 2 (p n/2 . (19) 

n=0 



In particular 



B ± = £ y 2 B i l /2) + e^ 2 B i l /2) + --- (20) 
B z = B + eB ( z l) + e 2 B {2) + ■■■ . (21) 



From the Faraday equation and the assumption of cold and massless electrons that implies E ■ B = 0, one has for the 
electric field 



E ± =e 5 / 2 E { l /2) + e y 2 E { l /2) + --- (22) 
E z = e 4 E^ ] + e 5 Ei 5) ---. (23) 

( 3 /2) 

One shows that the parallel current vanishes to leading order (Vj_ x B ± ' = 0) which, when combined with the 
divergenceless condition V ± •B ( | 3/2) + d z B z l) = 0, implies B ( | 3/2) = (-A ± )" 1 V ± 5 Z ^ 1) . 
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The ion velocity equation is then rewritten in the form of the pressure balance equation 

_ , B , , bl , 2e/ , m,. b z 



where, because of the cancelation of the leading order at threshold, we have defined b z — B z + and p± = 
p±^ +Ep^ .In this regime, the ion inertia term is subdominant, and the time dependency originates from the Landau 
resonance through the transverse pressure. This contrasts with the MHD approach where it arises from the advection 
term. 

In Eq. d24i >. the perpendicular pressure and the contribution of the gyroviscous tensor gradient are to be computed 

from the proton distribution function / = /(°) + e/' 1 ' + g 3 / 2 /( 3 / 2 ) + g 2 /( 2 ) -| , where the leading order is the 

equilibrium state and the other contributions are obtained by perturbative resolution of the Vlasov-Maxwell equations. 
Detailed computations are presented in [82] in the case of a bi-Maxwellian equilibrium state and in |)85T] for an arbitrary 
equilibrium distribution function. After returning to the unsealed variables, this leads to the dynamical equation 

(25) 



d t b = ^J |v(-jrd z ) (rb + ^r 2 A ± b-xd?Al l b-Ab 2 ^j , 



that involves several coefficients (and among them the distance to threshold T) whose expressions in terms of momenta 
of the equilibrium distribution and of its derivatives, are given in Ref. [85]- Here J$? denotes the Hilbert transform and 
the operator J(f z = —J^d z has the Fourier symbol \k z \. In the case of a bi-Maxwellian distribution function, Eq. d25l ) 



reduces to 

'2 



P±-P\\\>-^ , 3/1 



Equations ( 1251 ) or ( I26l l can be viewed as an extension of the linear dispersion relation by a nonlinear contribution that 
at the order of the expansion is not affected by FLR corrections, and can thus be estimated in the framework of the 
drift kinetic approximation [84, 86]. After a simple change of variables, Eq. d25l l becomes 

d z u = (-,^fd^)(au + A ± u-Aj 1 d^u-3sgn(A)u 2 ), (27) 

where a = sgnT is equal to 1 above threshold and —1 below threshold. Note that changing the sign of the nonlinear 
coupling coefficient A leads to replacing u by — u and thus magnetic holes by magnetic humps or the reverse. For a 
bi-Maxwellian equilibrium, A > 0. In this case, when looking for solutions depending on a unique space coordinate 
along a prescribed direction quasi-transverse to the ambient field, Eq. d26l i reduces to 



d t u = JSez[(ou + dz Z ) u-3u 2 )}. (28) 

A numerical integration of Eq. d28b above threshold with initial conditions in the form of a sinusoidal function of weak 
amplitude displaying several wavelengths in the integration domain, is presented in Ref. 18411 . After the early linear 
phase, several magnetic holes are formed, whose number is progressively reduced to one. Nevertheless, the asymptotic 
equation cannot capture the nonlinear saturation of the instability and, after a while, the solution blows up, indicating 
the breakdown of the asymptotic scaling and the formation of large-amplitude structures. Near collapse, the solution 
displays a self-similar profile discussed in lf84h . Furthermore, as detailed in Ref. l87ll . Eqs. d27i i and (E5T i do not admit 
stationary non-zero solutions above threshold, while below threshold steady solutions (in the form of magnetic holes) 
exist but are unstable. It indeed turns out that the instability is associated with a subcritical bifurcation, a property 
consistent with the bistability property reported from both satellite observations and numerical simulations. 



Nevertheless, the reductive perturbative expansion near a bi-Maxwellian equilibrium state, that predicts the forma- 
tion of large-amplitude magnetic holes above threshold, contradicts the results of direct numerical integrations of the 
Vlasov-Maxwell equations with similar initial conditions, which lead to magnetic humps. The question thus arises 
of the origin of this discrepancy which could be the consequence of one of the two main assumptions involved in 
the asymptotic theory, namely the fact that the ion distribution function remains close to a bi-Maxwellian, and also 
that kinetic effects are retained at the linear level only. In addition, as already mentioned, this theory does not capture 
the saturation of the instability which is in fact strengthened by the retained nonlinearity. Before revisiting the above 
theory, let us briefly review phenomenological models for the saturation of the mirror instability, previously discussed 
in the literature. 



6.2. Previous phenomenological modeling of the saturation 



The models proposed in [88, 89, l90ll are mostly based on a separation between trapped (with large pitch angle) and 
untrapped (with small pinch angle) particles which respond differently to magnetic field variations, the saturation of 
the instability resulting from the cooling of trapped particles in magnetic holes. In the rising field regions, trapped 
particles are indeed excluded by the mirror force, leading to a decrease of the particle pressure (reduction of j3jj and 
to the evolution to marginal stability (without important change in the particle energy). In the well regions, no particle 
can be excluded. Some trapped particles are cooled by loosing perpendicular energy, leading to a reduction of the 
temperature anisotropy. Large reductions in the field are required in the wells in order to cool the trapped population 
enough to stabilize the system. Such models mostly predict deep magnetic fields in conditions of marginal stability. 
They hardly explain the formation of magnetic humps (except possibly for unusually high /3 's), and do not address the 
phenomenon of bistability. 

Particle trapping is difficult to retain within a systematic reductive perturbative analysis. Within this formalism, it 
can nevertheless be phenomenologically described by a renormalization of the time derivative, on the basis that it 
produces a quenching of the Landau resonance. This argument is quantitatively implemented in ifs^tl by prescribing a 
flattening of the distribution function of the parallel velocity on a range that extends with the amplitude of the magnetic 
perturbation. This leads the authors to modify Eq. ( 1281 by multiplying the term d t u by a nonlinear function of u that 
prevents the explosive growth of the amplitude and, by vanishing, yields stationary solutions in the form of negative 
Korteweg-de Vries solitons. Yet, only magnetic holes can result from this model. 



6.3. Modeling the variation of the local ion Larmor radius 

As already mentioned, either the persistence to leading order of a bi-Maxwellian distribution function or the linear 
description of the kinetic effects (in particular the FLR corrections), that both result from the reductive perturbative 
expansion scalings, are possibly at the origin of the discrepancy with the results of the slightly supercritical direct 
numerical integrations of the Vlasov-Maxwell equations with an initial bi-Maxwellian ion distribution function. 

The model presented in l84ll consists in modifying Eq. d28l in order to phenomenologically take into account the 
variation of the local ion Larmor radius due to the development of magnetic structures. Indeed, in regions of weaker 
magnetic field, the ion Larmor radius is larger, making stabilizing effects of FLR corrections more efficient than in 
the linear regime. Consequently, the mirror instability is more easily quenched in magnetic field minima, making 
magnetic humps more likely to form in the saturating phase of the mirror instability. The conservation to leading order 
of the magnetic moment suggests that the square Larmor radius scales like the inverse longitudinal magnetic field, 
which leads to replace the FLR contribution Aj_u in Eq. d27i > by the nonlinear term rr^ Aj_« (possibly supplemented 
by a higher order correction in the form of a bi-Laplacian to improve accuracy). Here the parameter a is given by 



and is thus proportional to the (unsigned) bi-Maxwellian distance to threshold. Numerical 

integrations of the resulting equations are presented in l82ll . Figure 10 of this reference shows the evolution of the 
solution from an initial noise in the supercritical regime (a = 1.54). The singularity is arrested, the linear instability 
saturating with the formation of magnetic humps whose number decreases in time by an effect of coarsening that 
tends to strengthen the mirror structures. Very close to threshold (a = 0.05), the system evolves to a regime of ordered 
fluctuations with a slight dominance of the minimum amplitude. Below threshold (but not too far in the subcritical 
regime), large initial perturbations lead to the formation of magnetic holes. Figure 11, that displays the skewness of 



the solution versus aa (a measure of the signed distance to threshold weighted by a function of is qualitatively 
consistent with the observational data: the skewness is negative below threshold and becomes positive sightly above 
threshold. 



6.4. Modeling quasi-linear effects 

An important observation from the numerical simulations of the Vlasov-Maxwell equations near threshold concerns 
the flattening of the ion distribution function near the zero parallel velocity (Fig. 4 of M82I0 . that results from a 
diffusion in the velocity space possibly described by the quasi-linear theory. This approach (see e.g. Ir9lll . Section 
3.1) was in fact early applied to the mirror dynamics 111811 . It assumes space homogeneity and thus the absence 
of coherent structures. It can thus be valid at early times only. It requires also many modes in interaction, and 
consequently an extended domain. Furthermore, the extension of the quasi-linear theory to an aperiodic instability 
is not straightforward, but turns out to be justified in the case of the mirror instability near threshold because of 
its small growth rate <C ^v^m). The method can be sketched as follows. Let the instantaneous ion distribution 
function be written in the form f(x,v,t) — fo(v,t) + f\(x,v,t) where /o(v,f) = (f(x,v,t)) x corresponds to the space- 
averaged contribution. The perturbation f\ is decomposed in Fourier modes in the form f\ = \Y*k fk e>k x ■ Similarly, 
the electric field reads E = jLt^e'**. The Faraday equation prescribes the fluctuating Fourier modes in the form 
B/c = x E/c. Substituting in the Vlasov equation, one gets 

*A = -|;l{^( 1 -^) + 4 ( '' £;) } V - A+c " <29) 

A main hypothesis is that the variation of fo is slow compared with that of the fluctuations. Furthermore, the amplitude 
of the harmonic modes fa are assumed to be of small amplitude, and consequently their interactions are neglected. In 
this framework, the modes fa are expressed as linear functionals of fy. After some algebra, this leads to a diffusion 
equation in velocity space for the space-average distribution function that, when suppressing the subscript zero in order 
to simplify the writing, takes the form 

37 = 3 — ^llll^ ^ 3 — V -L [ D ±\\^— +Z) -L-L3 — (30) 
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with 
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supplemented by (bt — 8B z (Ic)/Bq) 

db k 

= K b k> ( 32 ) 

where Yk is the linear growth rate, computed using Eq. (O where the coefficients are evaluated using the instantaneous 
value of the space-averaged distribution function. The resulting system thus concentrates on the spatially independent 
part of the distribution function, ignoring the nonlinearities associated with space variations and describing the wave- 
wave interactions. It is solved numerically in II 8511 . The results are summarized in Fig. 1 of this reference. It clearly 
displays the saturation of the energy of the magnetic fluctuations. The system is shown to evolve towards a regime 
of marginal stability as the parameter T that measures the instantaneous distance to threshold decays to zero (an 
asymptotic value that is actually difficult to reach numerically). Simultaneously, the growth rate also decreases to 
zero. In addition, the gray scale plot of v±8f where 8f denotes the difference of the solution (corresponding to the 
space averaged distribution function) in the saturated regime with its initial value (Fig. 2) accurately reproduces the 
qualitative features of Fig. 4 of II 8 211 where the same quantity is evaluated from the direct numerical simulation of 
the Vlasov-Maxwell equations in an extended domain near threshold: both figures display a similar diffusion in the 
velocity space near the resonance. 



6.5. Coupling quasi-linear theory and reductive perturbative expansion 



The observation that the early nonlinear dynamics of the mirror modes is accurately described by the quasi-linear 
theory before the formation of strong coherent structures, suggests to couple the quasi-linear theory and the reductive 
perturbative expansion by estimating the coefficients entering Eq. d25l ) not from the (initial) equilibrium state but from 

the instantaneous quasi-linear distribution function ll85ll . a procedure consistent with the key assumption ^ -C Yk 
of the quasi-linear theory. Nevertheless, because of the nearly singular character of the distribution function resulting 
from the quasi-linear evolution near the zero parallel velocity, contributions of the resonant particles are to be taken into 
account when estimating the nonlinear coupling. This effects, which corresponds to the nonlinear Landau damping, 
leads to the denominator in Eq. (l34l 

df d df 1 d ( df df\ 

-37 = 3 — C|||| 3 1 3 — v ± [D ± n — +D ±± - — (33) 
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where the additional parameter vT 1 that measures the strength of the resonance, is given as an explicit functional of 
the instantaneous distribution function. The above system can also be viewed as retaining the full nonlinear evolution 
for the magnetic fluctuations given by the reductive perturbative theory, in the estimate of the coefficients of the quasi- 
linear diffusion equation, instead of the sole linear contribution. The resulting system was numerically integrated in 
the supercritical regime in ll85ll . One observes that while the time evolution of the distance to threshold and of the 
instantaneous growth rate are weakly perturbed compared with the genuine quasi-linear theory, the nonlinear coupling 
parameter A and the resonance coefficient v^ 1 change sign at early time (Fig. 3). Furthermore, after the linear phase, 
the local maximum amplitude of the magnetic fluctuations no longer saturate, but displays a fast growth associated 
with the formation of large-amplitude structures. As previously mentioned, the change of the sign of the A coefficient 
modifies the shape of the structures. Indeed, as displayed in Fig. 4, the emerging structures are now magnetic humps, 
as observed in the direct numerical simulations of the Vlasov-Maxwell equations. 

Capturing the saturation of the mirror structures by integration of Eqs. (T33l>— (f34b is rather delicate due to numerical 
limitations. In order to study the saturation effect produced by the denominator associated with the nonlinear Landau 
damping in Eq. (1341 . a simplified approach used in i92tl consists in freezing the coefficients after the early quasi-linear 
phase present in the supercritical regime, leaning on the quasi-singular character of the subsequent nonlinear dynamics. 
In the subcritical regime, the quasi-linear phase does not occur, but the denominator in Eq. fl34l l is nevertheless retained 
because of the large amplitude of the fluctuations needed to capture non trivial solutions. This leads to the model 
equation 

d,b = - ——^ (-Jfd^ )(ob + lidtfb- 3sb 2 ) , (35) 

where % is the coordinate along the most unstable direction. Here, a — ±1 depending on the super or sub-critical 
regime, and s = ±1 depending on the sign of the nonlinear coupling. The latter is negative above threshold after the 
early quasi-linear phase, and positive in the case of persistence of a bi-Maxwellian distribution as in the subcritical 
regime. 

Numerical integration of Eq. ( 1351 ) in the supercritical regime shows that during the nonlinear phase of amplitude 
growth, a plateau of negative values gradually develops, that tends to locally reduce the ambient magnetic field, putting 
the system in a situation similar to the subcritical regime. The solution is then attracted to a steady state given by 
the negative of the Korteveg-de Vries soliton profile and a maximal amplitude b max = 1/a. The amplitude of the 
structures is prescribed by the strength of the early-time quasi-linear resonance: larger amplitudes are obtained when 
these effects are smaller. When starting with random initial conditions, which leads to a large number of humps, a 
coarsening phenomenon is observed. If in the supercritical regime quasi-linear effects are subdominant, s — + 1 and 
magnetic holes are obtained. Below threshold (a — — 1), the distribution function remains bi-Maxwellian (s = +1) 
and large-amplitude magnetic holes are stable. 

The above discussion where the quasi-linear diffusion plays a main role, assumes an extended system. When because 
of geometrical constraints such as a small computational domain, only a few modes are unstable, particle trapping can 
play a role but its effect on on the distribution function is in fact similar to that of the quasi-linear diffusion, both 
processes leading to a flattening near the zero parallel velocity l86l l93ll and consequently to a change of the sign of 



the nonlinear coupling. A possibly more refined description of trapping could consists in computing the coefficients 
entering Eq. (f25l > on the basis of a distribution function that distinguishes between free ions, trapped ions and tail 
particles, as suggested in the construction of the steady mirror structures described in Ir94ll . 

7. TURBULENT PLASMA WITH MIRROR MODES 
7.1. Observational results 

Cluster measurements in the Earth's magnetosphere close to the magnetopause, presented in [95] concern a regime 
where the turbulence is compressible, anisotropic and dominated by mirror structures. Figure 1 of this letter displays 
high-level quasi-regular oscillations on the parallel component of the magnetic fluctuations, signature of strong 
compressible waves. The so-called k-filtering technique 119611 enabled the authors to discriminate between time and 
space variations. Figure 1 indeed shows a /~ 7 / 3 range for the low-frequency temporal spectrum, while Fig. 6 displays 
a k v 8//3 behavior for the spatial spectrum in the direction of the flow (referred to by the index v). It turns out that the 
"waves" are actually stationary structures in the flow frame, identified as mirror modes (see Fig. 3 and also I97l l4l0. It 
appears that energy is injected at large scales via a mirror instability well predicted by the linear theory, and cascades 
nonlinearly more preferentially along the flow direction (see Fig. 5), down to scales of the order of a few ion Larmor 
radii. 

The values of the proton temperature anisotropy and of the parameter of j3|| measured in the slow solar wind by the 
WIND/SWE mission are considered in [98], and their relation with the temperature-anisotropy-driven instabilities are 
discussed. A main result is that the proton temperature anisotropy is constrained by oblique instabilities, either mirror 
or firehose, depending on the dominance of the perpendicular or parallel temperature respectively. This observation 
contrasts with the linear theory which predicts a dominance of the proton cyclotron and of the parallel fire-hose 
instabilities. It is exemplified in Fig. 1 of ll98ll that shows a color-scale plot of the relative frequency of (fi\\ p , Tj_ p /Ty) 
for the considered data, together with contours of estimated growth rates of the above instabilities in a plasma with 
bi-Maxwellian ions and isothermal electrons. A similar conclusion is reached in B99I1 by analyzing 10 6 independent 
measurements of gyroscale magnetic fluctuations in the solar wind. An additional result is that these fluctuations are 
enhanced along the temperature threshold of the mirror, proton oblique firehose, and ion cyclotron instabilities. In 
addition, the relative importance of the parallel magnetic fluctuations, as measured by \8b\\ \ 2 / (\8b« \ 2 + \8bj_\ 2 ), is 
increased at j8|| > 1 along the mirror threshold but small elsewhere, as expected for a mirror instability. 

The question then arises of the origin of the temperature anisotropy. As already mentioned, in the planet mag- 
netosheaths large-perpendicular-to-parallel ion temperature anisotropy ratio are due to compression downstream the 
quasi-perpendicular shock or to the flow of plasma with high parallel velocity close to the magnetopause. In the solar 
wind, the mechanism that is usually believed to drive such an anisotropy is the proton resonance with left-handed 
cyclotron waves. However, as recently noticed in |6], this process requires a parallel cascade via parametric decay to 
high frequency modes, while the development of quasi-perpendicular low-frequency kinetic Alfven waves is usually 
considered to be the dominant process. The statistical analysis presented in 16J] of data in the fast solar wind provided 
by Helios 2 primary mission at radial distance from the sun from 0.3 to 0.9 AU (where J3|| ranges between about 0.1 
and 0.6), shows a correlation between enhancements of the power density of transverse waves (with frequency be- 
tween 0.01 and 0.1 ion gyrofrequency in the plasma frame) and increases of the ratio of the perpendicular-to-parallel 
proton temperatures. This led the authors to conclude that perpendicular ion heating is connected with the presence of 
a large-scale turbulence where transverse oscillations are dominant. 



7.2. Numerical simulation of the mode competition 

It has often be argued 1T251 l69ll that, in a plasma with a dominant perpendicular temperature, the presence of He ++ 
which lower the linear ion-cyclotron growth rate without significantly affecting the mirror modes, could be at the origin 
of the dominance of the latter structures in planetary magnetosheaths. Such an ambiguity associated with multiple ion 
effects is not present in numerical simulations where an electron-proton plasma can be considered. Such simulations 
were performed in J5t] using hybrid PIC simulations in two and three space dimensions. Figure 2 of this paper shows 
that in the linear regime (see also [30), the ion-cyclotron (referred to as L-mode EMIC) growth rate dominates that 
of the mirror modes. This ordering does not necessarily persist in the non linear regime. To address this question, the 



authors identify the magnetic energy density of each wave by the following method. Examination of the wave spectra 
shows that ion cyclotron waves and mirror mode waves are separated by the angle 9 = 30° between the wave vector 
k and Bq. Accordingly, the "L-mode EMIC wave range" is assumes to be defined by 0° ^ 9 ^ 30°, and the "mirror 
mode wave range" by 30° < 9 ^ 80°. The spectra of the magnetic field is divided into two ranges of the wave vector 
space, and the energy associated with the two types of waves is obtained by integrating the square magnetic modes 
|B(fe)| 2 /|Z?o| 2 in each of these ranges. Figure 6 compares the time evolution of the energy associated with the two 
types of waves, both in two and three dimensions. It appears that in two dimensions, the mirror energy is strongly 
subdominant. Differently, in three dimensions the two energies reach the same maximum after about the same time 
but, while the mirror energy then remains quasi-stationary, the ion-cyclotron energy rapidly decays in time making the 
mirror modes strongly dominant in the long-time regime. 



7.3. Toward fluid simulations of the mirror dynamics 

Strictly speaking an accurate description of collisionless plasmas even in the presence of a strong ambient field 
requires a fully kinetic description, but in the turbulent regime huge computational resources are requested for such 
simulations. On the other hand, it is clear that usual MHD or bi-fluid descriptions are unable to accurately describe such 
a regime, as exemplified by the mirror instability. It is thus of interest to extend the fluid approach in order to include 
leadi ng low-frequency kinetic effects such as Landau damping and FLR corrections. This approach was initiated in 
[100] where a closure relation retaining linear Landau damping was suggested. It was implemented in the case of 
large-scale MHD in [77], as a closure of the moment hierarchy in the framework of the kinetic drift approximation, 
leading to a s o-called L andau fluid model. It was extended to dispersive MHD with Hall effect and large-scale FLR 
corrections in generalization of the model including a description of quasi-transverse scales significantly 

smaller than the ion gyroscale, under the gyrokinetic scaling, is presented in 112611 and referred to as FLR-Landau 
fluid model. Note that in contrast with the gyrokinetic description, Landau fluids retain fast waves that are accurately 
described up to the ion gyroscale. Nevertheless, like the gyrofluids, they neglect wave-particle trapping, that is to say 
the effect of particle bounce motion on the distribution function near resonance. 

We cannot detail here the derivation of the FLR-Landau fluid model nor even write the full equations explicitly. In 
brief, it consists in deriving from the Vlasov-Maxwell equations for the protons and the electrons a hierarchy of exact 
equations for the plasma density and velocity (for the sake of simplicity the electron inertia can be neglected, leading 
to a generalized Ohm's law), and for the gyrotropic pressures and heat fluxes of each particle species. The closure 
concerns the gyrotropic components of the fourth order cumulants of each species and the non gyrotropic components 
(FLR corrections) of the pressures, heat fluxes and fourth order moments of the ions (electron FLR corrections are 
neglected, as the description does not go beyond the ionic scales). It is built in such a way that the linear kinetic theory 
is reproduced in the low-frequency limit, within the gyrokinetic scaling in the quasi-transverse directions but also at 
all the hydrodynamic scales. For this purpose, the various moments are first computed within the usual linear kinetic 
theory in the above asymptotics, leading to expressions in terms of the electromagnetic fluctuations. The difficulty is 
that these formulae involve the electrostatic response function (or the plasma dispersion function) that depends on the 
phase velocity of the corresponding Fourier mode and consequently cannot be implemented in the framework of an 
initial value problem. The idea is first to combine these expressions for the various moments in order to express the 
unknown cumulants in terms of lower order moments and of the magnetic current. In such expressions, the response 
function arises within a rational function. At this step, the response function is replaced by a Pade approximant whose 
number of poles is chosen in such a way that, when returning to the physical space, the closure relations reduce to 
partial differential equations, usually of the first order in time. They involve an Hilbert transform in the longitudinal 
space variable, signature of the Landau resonance. Being derived in the context of the linear kinetic theory, these 
closures involve equilibrium quantities such as the parallel and perpendicular pressures or temperatures. In order to 
take into account the global evolution of the plasma, such as its global heating, it is suitable to replace these equilibrium 
(i.e., initial) quantities by the instantaneous space-averaged values. A similar approach is implemented to evaluate 
the relevant FLR corrections. A detailed derivation of the model is presented in [26]. Explicit expressions for the 
various contributions in the one-dimensional regime where the space variations only take place in a direction making 
a prescribed angle with the ambient field, are given in ill 0311 . 

Detailed comparison with the kinetic theory demonstrates that the FLR-Landau fluid model correctly reproduces 
the mirror instability and its quenching at small scale (Fig. 10 of [26]). It also accurately reproduces the frequency and 
the damping rate of the (quasi-transverse) kinetic Alfven waves (Fig. 2). In oblique directions (Fig. 3), the accuracy is 



limited by the presence of resonances. 

The nonlinear dynamics of mirror modes in one space dimension is studied using the FLR-Landau fluid in Il03ll . 
Compared with Eq. (l28l which results from the reductive perturbative expansion of the Vlasov-Maxwell equations, 
the FLR-landau fluid model involves richer nonlinearities. As a consequence, it is capable to arrest the singularity. 
Nevertheless, in both models, the distribution function remains close to a bi-Maxwellian and FLR corrections are 
described at an essentially linear level. As a consequence, in both cases, the saturation of the mirror instability leads to 
magnetic holes and not humps. Figure 2 of lll03ll exemplifies in the FLR-Landau fluid framework the development of 
the mirror instability and the formation of magnetic holes antic orrelated with density humps, and their slow relaxation 
towards the uniform state on a long time scale of the order of 10 fir 1 , with the temperatures evolving toward a 
mirror stable regime (Fig. 3). In order to prevent this relaxation, simulations were also performed when the mean 
ion pressures are prescribed to constant values. In a computational domain of relatively small extension, the system 
evolves to a stationary state (Fig. 4). Through the magnetic holes, the increase of the perpendicular temperature is 
larger than that of the parallel one. It is of interest to note the agreement with observations reported in [ 104] that solar 
wind protons detected within magnetic holes are found to be preferentially heated perpendicular to the ambient field. 
Simulations of the FLR-Landau fluid with prescribed mean pressures in an extended domain and at moderate distance 
from threshold (Fig. 5) display a regime of spatio-temporal chaos. Holes form and disappear in an erratic fashion, and 
large-scale compression waves are observed to propagate through the domain. 

Furthermore, below the mirror instability threshold, large-amplitude stationary solutions have been obtained using a 
continuation procedure (Fig. 6). Moreover, the solution displays a small magnetic field component perpendicular to the 
(fc,2?o)-plane (Figs. 7 and 8). Interestingly, this component is symmetric with respect to the center of the magnetic hole, 
a property also obser ved in direct Vl a sov- Maxwell simulations ((HI, Fig. 13), which contrasts with soliton models 
based on Hall-MHD lll05l Il06l Il07l Il08ll . All the velocity components are anti-symmetric, while the component 
perpendicular to the (k,Bo)-plane has a relatively large amplitude with a sharp gradient at the center of the hole (Fig. 8), 
indicating that the gyroviscous tensor plays an important role in the equilibrium. Similar signatur es are observed in 
PIC simulations of non-propagating rarefractive solitary structures generated by particle injection II 1 0911 in a plasma 
that is initially isotropic. The heat fluxes also seem to play an important role since no equilibrium can be reached when 
these terms are removed from the equations. 

Although the FLR-Landau fluid model does not accurately reproduce the geometry of the nonlinear mirror modes in 
the weakly unstable regime, the influence of the mirror instability on the development of the temperature anisotropy of 
a turbulent plasma can nevertheless be correctly captured. In this context, it is of interest to mention recent simulations 
performed in one space dimension along a direction quasi-transverse to the ambient field 111 1011 . A turbulent regime is 
considered where the system is driven by a random force acting when the sum of the kinetic and magnetic energies 
is smaller than a prescribed value (changing this value mostly affects the evolution time scale), in a way aimed to 
create large-scale transverse oscillations. It turns out that for J3|| larger than a critical value close to 0.35, a temperature 
anisotropy corresponding to a dominant perpendicular ion heating develops. As time elapses, the system reaches a 
regime where the mirror instability develops and constrains the anisotropy at a level slightly in excess of the threshold 
condition estimated by Eq. ©. In the physical space, mirror structures are clearly identified. At smaller j3||, the parallel 
ion temperature dominate, leading to significantly different dynamics. More detailed analyses of these runs are in 
progress, but the present simulations already provide an interesting insight on the role of turbulence a nd non-res onant 
interactions on the temperature anisotropy observed in the solar wind, whose origin is still debated 13, 1 1 1 ll[lT2tl . 

8. SUMMARY 

We have presented an overview of the present understanding of the dynamics of the mirror modes commonly observed 
in space plasmas displaying a dominant ion perpendicular temperature and a relatively large j3 . Various questions arise, 
originating in particular from the existence of mirror structures in plasmas that are stable for the mirror instability, and 
also from their appearance as both magnetic enhancements or depressions, with a statistical predominance of the latter. 
A main conclusion resulting from recent analyzes of the mirror data provided by satellite missions is that magnetic 
humps are only observed in an unstable plasma, while in the mirror-stable regime only magnetic holes are present. 
Direct numerical simulations of the Vlasov-Maxwell equations for initially bi-Maxwellian ions (and cold electrons) 
display a flattening of the ion distribution function near the zero parallel velocity, viewed as the result of a quasi-linear 
diffusion in the velocity space, together with the development of persistent magnetic humps. At larger distance from 
threshold and for a j8|| that is not too large, the flattening of the distribution function is not observed and, after a while, 
the early formed magnetic humps can evolve toward magnetic holes. A similar change of shape is obtained when 



retaining the expansion of the plasma, which makes the system to progressively evolve to a stable regime. Below 
threshold, magnetic humps are damped but magnetic holes turn out to be persistent. This suggests that the magnetic 
holes displaying the polarization properties of the mirror modes, result from the evolution of magnetic dips and humps 
possibly convected to the point of observation. 

A fully systematic theory reproducing all the observational and numerical results is still to be developed. A reduc- 
tive perturbative expansion leads to an asymptotic equation that can be viewed as an extension of the growth rate 
equation, retaining the leading nonlinearity with a coupling coefficient whose sign prescribes the shape (depressions 
or enhancements) of the mirror structures. In the case of a bi-Maxwellian initial equilibrium, this approach predicts 
magnetic holes. It turns out that taking into account the early quasi-linear evolution changes the sign of the nonlinear 
coupling coefficient, and thus reproduces the formation of the magnetic humps observed in the direct numerical simu- 
lations. The reductive perturbative expansion also predicts that the bifurcation associated with the mirror instability is 
subcritical, which provides an interpretation of both the phenomenon of bistability (i.e., the stability of mirror struc- 
tures below threshold) and the large amplitude of the observed mirror structures. The latter property is associated with 
the development of a finite-time singularity for asymptotic solutions above threshold, indicating a breakdown of the 
asymptotic scaling. Nonlinear Landau damping provides a plausible mechanism to saturate the structure amplitude. 

In addition to their fundamental interest, linear and nonlinear mirror modes are expected to affect the transport 
properties of the plasma and also to constrain the temperature anisotropy associated with a dominant perpendicular 
heating of the ions, observed in the solar wind. As recently suggested by simulations of a refined Landau fluid model, 
such a heating could result in particular from the quasi-transverse dynamics of a turbulent plasma driven by large-scale 
quasi-incompressible oscillations, a scenario that is consistent with the classical picture of a quasi-transverse Alfvenic 
cascade and with previous suggestions of a non resonant ion heating by Alfven wave turbulence. 
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